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Two geometrical well-posed hyperbolic formulations of general relativity are de- 
scribed. One admits any time-slicing which preserves a generalized harmonic con- 
dition. The other admits arbitrary time-slicings. Both systems have only the 
physical characteristic speeds of zero and the speed of light. 

Einstein's theory, viewed mathematically as a system of second-order par- 
tial differential equations for the metric, is not a hyperbolic system without 
modification and is not manifestly well-posed, though physical information 
propagates at the speed of light. A well-posed hyperbolic system admits 
unique solutions which depend continuously on the initial data and seems 
to be required for robust, stable numerical integration. The well-known tradi- 
tional approach achieves hyperbolicity through special coordinate choices. The 
formulation£l described here permit much greater coordinate gauge freedom. 
Because these exact nonlinear theories incorporate the constraints, they are 
natural starting points for developing gauge-invariant perturbation theory. 

Consider a globally hyperbolic manifold of topology Exi? with the metric 

ds 2 = -N 2 di 2 + gi j (dx i + f3 l dt)(dx 3 + f3 j dt), (1) 

where N is the lapse, (3 Z is the shift, and is the spatial metric. The deriva- 
tive, N do = N~ 1 (d/dt — Cp), where Cp is the Lie derivative in a time slice 
X along the shift vector, is the derivative with respect to proper time along the 
normal to S. This implies that do is the natural time derivative for evolution. 
The extrinsic curvature Kij of S is defined as 

d ogij = -2NK ir (2) 

Together with (|^), the dynamical part of Einstein's theory (in vacuum) 
can be expressed in 3+1 language as 

Rij = -N^Ky + + Rij = 0, (3) 
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where Rij is the spatial Ricci curvature and jij consists of terms at most zeroth 
order in derivatives of K, first order in derivatives of g, and second order in 
derivatives of N. The spatial Ricci curvature is second order in derivatives of 
g in such a way as to spoil the hyperbolicity of the combined equations (|J) and 
(||). To achieve hyperbolicity for the 3+1 equations, we proceed as follows. 

By taking a time derivative of (JsJ) and subtracting appropriate spatial 
covariant derivatives of the momentum constraints, we obtain an equation 
with a wave operator acting on the extrinsic curvature. In vacuum, one finds 

doRij - ViRoj - VjRoi = NnKij + J l3 + % = 0, (4) 

where □ = —N~ 1 doN~ 1 do + V fc Vfe, Jij consists of terms at most first order 
in derivatives of K, second order in derivatives of g, and second order in 
derivatives of N, and 

Sij = -N-^iV^doN + N 2 H) (5) 

(H = K k). The term Sij is second order in derivatives of K and would spoil 
hyperbolicity of the wave operator □ acting on K. Hyperbolicity is achieved 
by imposing "generalized harmonic slicing" 

d a N + N 2 H = g^doig-^N) = f(x, t; g- 1/2 aN, g 1 ' 2 ^ 1 ), (6) 

where / is an arbitrary scalar function of its arguments and a is an arbitrary 
scalar density of weight one. 

After using (^) to replace K in (H) and designating gij and g~ x l 2 uN as 
variables, the resulting equation and (0) form a quasi-diagonal hyperbolic sys- 
tem with principal operators <9 £l and dg. This system can also be put in first 
order symmetric hyperbolic form by introducing sufficient auxiliary variables 
and by using the equation for i?oo (thus incorporating the Hamiltonian con- 
straint). The Cauchy data for the system (in vacuum) are: i) (g, K) such that 
the constraints i?oi = 0, G°q = hold on the initial slice; ii) 9oK such that 
Rij — on the initial slice; and iii) N > arbitrary on the initial slice. Note 
that the shift f3 k (x,t) is arbitrary. Using the Bianchi identities, one can prove 
that this system is fully equivalent to the Einstein equations. The constraints 
of this system propagate in a well-posed causal hyperbolic way. The restric- 
tion on the initial value of <9oK is what prevents the higher derivative from 
introducing spurious unphysical solutions. 

All variables propagate either with characteristic speed zero or the speed 
of light. The only variables which propagate at the speed of light have the 
dimensions of curvature, and one sees that this is a theory of propagating 
curvature. 
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In the above formulation, the shift is arbitrary, but the lapse is determined 
by (^|). The theory is therefore covariant under coordinate transformations 
which preserve the slicing. Can the lapse be made arbitrary so that one can 
simply drop the generalized harmonic condition (^)? The answer is yes. 

By taking another time derivative and adding an appropriate derivative of 
i?oo, one finds (in vacuum) 

d d Rij - doViRoj - ^Vj-JZw + %VjRoo = MNQK^) + J y = 0, (7) 

where Jij consists of terms at most third order in derivatives of g and second 
order in derivatives of K. Together with (|J), these forrrLa system for (g,K) 
which is hyperbolic non-strict in the sense of Leray-OhyaB Both the lapse and 
the shift are arbitrary (N > 0) and are not dynamical variables. The Cauchy 
data of the previous form (in vacuum) must be supplemented for (||) and 
by (9o<9oK such that d$Rij = on the initial slice. This guarantees that the 
system is fully equivalent to Einstein's theory. This system does not have a 
first order symmetric hyperbolic formulation. 
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